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m ; Abstract 

o 

' The method of the effective action for the composite operators <I>^(x) 

and is apphed to the termodynamics of the scalar quantum field 

Q\ i with interaction. An expansion of the finite temperature effective 

potential in powers of h provides successive approximations to the free 
energy with an effective mass and an effective coupling determined by 
C^. the gap equations. The numerical results are studied in the space- 

(— I I time of one dimension, when the theory is equivalent to the quantum 

mechanics of an anharmonic oscillator. The approximations to the 
free energy show quick convergence to the exact result. 
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1 Introduction 



Thermal properties of quantum field systems are important for a study of 
phase transitions in the big bang cosmology and for an interpretation of heavy 
ion collisions. These systems are characterised by a large coupling constant, 
thus non-perturbative methods should be developped to study their termo- 
dynamics. Here we shall discuss an application of the effective action for 
local composite operators to the finite temperaturequantum field theory. 
We consider a scalar field with a classical action in n- dimensional Euclidean 
space-time given by 

Sm = j^^x){-d^ + m^)^{x) + \^\x)] rx, (1) 

but the method can be easily extended to fermionic systems. 

Recently, the diagrammatic rules for the effective action for the operator 
have been found and the method has been extended to include the 
operator $^(2;) 0. The lowest order approximation coincides with the Gaus- 
sian Effective Action resulting from the variational principle with Gaussian 
trial states [Q. Within the thermal field theory formalism this would corre- 
spond to the finite temperature Hartree approximation for the free energy 
An expansion of the effective action for local composite operators in powers 
of h provides thus a non-perturbative method for a systematic improvement 
of the Hartree approximation both at T = and T ^ 0. An advantage of 
using local operators lies in the fact that calculation is simpler than for bilo- 
cal composite operators ^{x)^{y) and other systematic methods which 
contain the Gaussian Effective Action as the lowest approximation (f.e. the 
optimized expansion 0). 

2 Finite temperature effective potential for 
composite operators 

The vacuum functional for the composite operators ^"^{x) and may be 

represented by a path integral 

Z[J, K] = = / D$ eM-^[*]+5 fKix)^^^) d-x] 
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and the effective action is defined through a Legendre transform 



r[A,A] = W[J,K]-^ lj{x)A{x)rx-^jK{x)A{x)rx 
r 

- — K{x)A\x)rx, (3) 
8 J 

where the background fields are given by 

fi'A(i) = 24-^5--3fi=A=(i)=<4.V)>J,*--3<<l''W>5K (4) 

oKyxj 

and < ... >j^K denotes the expectation value in the presence of external 
sources J and K. The physical quantities have to be calculated at J = K = 0, 
or equivalently at the values of the background fields for which the gap equa- 
tions 

6r h ^, ^ 



5A{x) 2 ^ ' 4 



-J{x) - —A{x)K{x) = 0, 

K{x) = (5) 



5A{x) 24 
are satisfied. 

The path-integral quantization in Euclidean space-time enables us to 
study the quantum field theory and its equilibrium thermodynamics by the 
same "imaginary time" formalism, if the appropriate boundary conditions are 
chosen ^. The vacuum functional (Q) defines the generating functional 
for composite Green's functions in quantum field theory, if the path integral 
is taken over the functions which approaches some constant as Xq ±oo. If 
the set of functions is periodic with a period = |t in the "imaginary time 
coordinate" xq, the integral @ has a meaning of the partition function at 
the temperature T, generalised to include interactions with the external cur- 
rents J and K. The finite temperature effective potential for the composite 
operators $^(2;) and is defined as 



(6) 

A.=const,A=const 
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The free energy density can be obtained as Fp = Vp{Aj3,A[^), where the 
temperature dependent expectation values of the composite fields, and 
A^, which correspond to J = K = 0, are determined by a stationary point 
of ^^(A, A), according to (||). 

Representing the connected generating functional by the series 

oo 

W[J,K] = Y.h'W^k)[J,Kl (7) 

k=0 

the effective action for composite operators (|^) may be obtained by eliminat- 
ing the sources 

oo oo 

J=^nV(fc)[A,A], and K = ^ n'=K(fe)[A, A] (8) 

k=0 fc=0 

in favour of the expectation values of the composite fields 

oo oo 

A = ^;i'=A(fc)[J,K], and A = ^n%)[J,K], (9) 

k=0 k=0 

order by order in h. The coefficients Jq and Kq are determined implicitly by 
the lowest-order relations 




A{x) = Ao{x) = G{x,x) = (x) (b j ^]_0) 
and 

A{x)=Aoix)=JG\x,y)iKoiy)~2AX)ry, (11) 

where the inverse propagator is given by 

G-\x,y) = {-d' + Q\x))6{x-y), (12) 
with an effective mass fl{x) defined by 

n^[A] = - Jo[A]. (13) 
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Higher coefficients, J(fc) and K(^k) niay be easily expressed as functionals of 
Jo [A] and iCo[A,A] and the Legendre transform @ can be performed. To 
the order h'^ the result is 

r[A, A] = J[^{n\x) - m2) A(x) - 3Xh^A\x) - h^XA{x)]dx 
with Q related to A by ([ToD and (0), and the effective coupling 




where a slash denotes an inversion of the operator. 

The effective action which includes the two-particle operator only, 
r[A], can be obtained from the above expression for r[A, A] by setting 1^0 = 
which results in 

r[A] = - m2)A(x) - 3\h^A\x)]dx 

in agreement with diagrammatic rules proven in Ref. [@. 

The expansions of r[A] and r[A,A] in powers of h provide two differ- 
ent approximation schemes for physical quantities, we shall denote these as 
12 and 14, respectively. The non-perturbative character of the approxima- 
tions is due to the fact that the expectation values of the composite fields are 
determined from the gap equations, obtained by requiring the given order ex- 
pression for the effective action to be stationary. At order both approaches 
results in the Gaussian effective action which contains an infinite sum of 
bubble diagrams of perturbation theory. As it happens, the approximations 
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to physical quantities obtained from r[A] and r[A, A] coincide also at order 
h^, after solving the gap equations, and differences first appear at order h'^. 
One can note that by introducing higher composite operators new approxi- 
mation schemes can be formulated, but their influence will appear in much 
higher orders of h. 

The finite temperature effective potential for composite operators can be 
obtained from the effective action by replacing the Feynman rules by those at 
finite temperature. In the given order approximation to the effective action 
the infinite interval in all integrals over t is replaced by [0,/5]. A and A 
are taken to be space-time independent, in this case Q is also constant and 
the Fourier transform of the propagator (0) can be performed. Using a 
momentum representation in space, but a coordinate representation in time 
we have 

where Um = are the Matsubara frequencies. The free energy calculated 

at the order h"^, appears identical with that in the Hartree approximation, 
provided by the finite temperature Gaussian Effective Potential The 
post-Gaussian corrections can be determined by taking into account higher 
orders of the finite temperature effective potential for composite operators. 



3 Free energy of the quantum-mechanical an- 
harmonic oscillator 

In the space-time of one dimension (time) the theory is equivalent to a 
quantum mechanical anharmonic oscillator with a Hamiltonian 

H = ^p'^ + ^m^x^ + Xx^. (18) 
After rescaling all quantities in terms of A, only one dimensionless parameter 

2 

z = ^2/3 remains; therefore, when discussing numerical results, we put A = 1 
without a loss of generality. The spectrum of the anharmonic oscillator can 
be calculated numerically and provides the simplest test for approximation 
methods in quantum field theory. The free energy, which contains informa- 
tion on the whole spectrum, can be used to study the reliability of field- 
theoretical methods at finite temperature. The conventional loop expansion 
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for the free energy coincides with termodynamic perturbation theory [p!0|] , 
the non-perturbative methods provide better approximations. 

Approximations to the ground state energy and to the second and fourth 
excitation of the anharmonic oscillator, calculated with the use of operators 

and have been shown to converge quickly to the exact results 

we shall discuss therefore the approximations to the free energy calculated 
from the termal effective potential for these operators. The temperature 
dependent propagator ([T7|) in one dimensional space-time becomes equal to 

' fisinh[f ] ^ ' 

and the integrals over /3 > t > in Feynman diagrams can be easily per- 
formed. We shall compare two approximation methods: 
(12) using the operator $^(a;): 

Fp = V^(A^), where V^(A) is calculated from Eq. ([T6|) and is 
a solution of the gap equation 



SA{x) 

(14) using the operators ^^{x) and $^(x) 



(20) 



Fjs = Va(A^, A^), where Va(A, A) is calculated from |I^ and A/3 
and Af3 are obtained as a solution of the gap equations 



SVr- 



6A{x) 



and 



5Vf 



5K{x) 



0. 



(21) 



The non-perturbative character of the approximations to the free energy 
is kept by solving the gap equations non-perturbatively (if one expanded 
the solutions of the gap equations to the given order in Ti, the results would 
coincide to that order in A with that obtained in thermodynamic perturbation 
theory). In both cases we choose the solution of the gap equations with the 
largest positive value for VL. The free energy for three values of the parameter 

2 

z = ^j3, z = 5,0,-1, is shown as a function of the inverse temperature 
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P = ^ in Figs. (l)-(3), respectively. The results obtained from Va(A) in 
tree lowest orders of h (121, 122 and 123) are compared with that obtained 
from Vg(A,A) at order (143). The exact free energy was obtained as 
= lnJ2n e~^^" where the energy levels, have been calculated by a 
linear variational method, using harmonic oscillator wave functions with an 
appropriately chosen frequency [0. 



The comparison shows that an expansion of the termal effective potential 
for composite operators in powers of h provides a method to improve the fi- 
nite temperature Hartree approximation (121) in a systematic way. Both the 
methods, 12 and 14, work very well for the single well anharmonic oscillator, 
the last being better in the whole range of temperatures. In both cases a 
quality of approximations becomes worse for decreasing values of z and in 
the double well case {z < 0) the methods are applicable only if the wells 
are not too deep {z > — 1). For z = —1 (Figure 3) large discrepancies 
between the results of different orders of the method 12 and the exact value 
appear. The results obtained at order with the use of the operators $^(x) 
and $^(a;) (143) are definitely better than that obtained with the operator 
$^(x) only (123). At low temperature discrepancies do not exceed 2%, but 
they grow with increasing temperature, where a contibution of higher excita- 
tions becomes more important. This can be understood from the analysis of 
the method 14 for the spectrum of the anharmonic oscillator , where it has 
been demonstated that the quality of the approximations becomes worse and 
worse for increasing excitation level. There is a hope that including higher 
composite operators will broaden the range of a good convergence further, 
making the method applicable even in the case of the double well potential 
with deeper wells. It will be also interesting to discuss an application of the 
local composite operators method to study thermodynamics of the scalar 
quantum field theory in the case of untrivial space dimension, we reserve this 
matter for a future publication. 
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Figure captions 



Figure 1. The free energy of the anharmonic oscillator ai z = ^ = 5, plot- 
ted vs p ^ ^. The successive orders approximations obtained from V^[A] 
{dashed line, 121, 122, 123) and from Vg[A,A] (dotted line, 143) compared 
with the exact value [solid line). 

Figure 2. As in Fig.l, but at z = 0. 

Figure 3. As in Fig.l, but at 2; = —1. 
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Figure z 
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